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Abstract 

The continuity equation relating the change in time of the position probabihty den- 
sity to the gradient of the probabihty current density is generahzed to PT-symmetric 
quantum mechanics. The normahzation condition of eigenfunctions is modified in ac- 
cordance with this new conservation law and illustrated with some detailed examples. 
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1 Introduction 



From time to time, non-Hermitian Hamiltonians have found applications in several research 
areas, such as nuclear physics, quantum field theory, condensed matter physics, and biol- 
ogy. Since the pioneering works of Bessis |l| and Bender and Boettcher 0, a subclass of 
such Hamiltonians, containing operators invariant under joint action of parity (P: x —x, 
p —p) and time reversal (T: x ^ x, p ^ —p, i — > — i), has become a subject matter 
of considerable interest. One of the main reasons for this is that there is a strong ana- 
ljd;ical and numerical evidence supporting the conjecture that, except when PT symmetry 
is spontaneously broken, the bound-state eigenvalues of these Hamiltonians are real (see 
e.g. i-ll). 

PT-symmetric quantum mechanics is often seen as a testing bench for nonconventional 
quantum field theories wherein PCT invariance (or PT invariance in scalar theories) plays a 
leading role (see e.g. |jT3|, [l^). This may be compared to the introduction of supersymmetric 



quantum mechanics by Witten as a testing ground for non-perturbative methods of 
achieving supersymmetry breaking in field theory. 

Both the field-theoretic interpretation and the phenomenological relevance of complex 
potentials in quantum mechanics make it necessary to better understand some fundamental 
issues connected with the replacement of Hermiticity by PT-symmetry, such as the physical 



interpretation of the Hamiltonian eigenfunctions |]T6|, [T^, |20|| . 

The present work is intended as a contribution to such an understanding. Our purpose 

is twofold. Firstly, we will generalize to PT-symmetric quantum mechanics the conservation 

law of standard quantum mechanics (here formulated for one-dimensional systems) 

dP{x,t) dJ{x,t) _ 

dt ^ dx ' ^ ' 

similar to the continuity equation of classical hydrodynamics and relating the change in 

time of the position probability density P{x, t) = \tp{x, t) p to the gradient of the probability 

current density J{x,t) = (fi/2mi) [ip*{x,t){dilj{x,t)/dx) — ip{x,t){dilj*{x,t)/dx)]. Secondly, 

we will discuss with some detailed examples a proposal for modifying the definition of 

the normalization condition of bound-state eigenfunctions in accordance with this new 

conservation law. 
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2 Generalized Continuity Equation and Some of Its 
Consequences 



Let us start with the time-dependent Schrodinger equation for a single particle in one 
dimension, 

(2) 



ifi — = — tt^ h V[x)ip[x,t), 



dt 2m dx"^ 

where the potential V{x) is assumed to be complex and PT-symmetric, i.e., PTV(x)PT 
V*{—x) = V{x). As a consequence, the function ip*{—x,t) satisfies the equation 

dij*{-x,t) fi2 d2^*(^-x,t) 



- ih- 



+ V{x)^lj*{-x,t), 



(3) 



dt 2m dx'^ 

where the potential is the same as in (H). 

Consider now eq. (0) for some solution ipi{x,t) and eq. (H) for some other solution 
'ip2{x, t), which may be equal to or different from ipi{x, t). Multiplying (|2|) by ■?/'2(— x, t), (0) 
by ilJi{x,t), and subtracting, we obtain 



|[^*(-x,tM(x,t)] + ^|- 



0. (4) 



dx ' dx 

For ip2{x,t) = ipi^x.t) = ip{x,t), it is clear that eq. (H) is a natural generalization of the 
continuity equation (||) of standard quantum mechanics, namely 

dPpT{x,t) , dJpT{x,t) 



dt 



dx 



0, 



(5) 



with 



Ppt{x, t) = ij*{-x, t)ilj{x, t), Jpt{x, t) 



2mi 



dx 



dx 



6) 



If ilJi{x,t) —>■ and ^/'2(x, t) ^ for x ^ ±oo, as is normally required for bound-state 
wave functions, integration of (HI) over the entire real line leads to the conservation law 



d r+co 
di 



dx ip2{~x,t)ipi{x,t) = 0, 



(7) 



or for 4'2{x, t) = ^ipii^x, t) = ip{x, t), 



— / dx PpxixA) = 0. 



Note that for Hermitian Hamiltonians, the result corresponding to is interpreted in 
standard quantum mechanics as the conservation of norm in time pT| . 

From now on, we are going to restrict ourselves to energy eigenfunctions 

t) = ui{x)e-^^'\ ij2{x, t) = U2{x)e~^^^\ (9) 

associated with (possibly complex) energies Ei, E2, respectively. On inserting (H) into (|^), 
we obtain 

(^1 - E;)e—n^^'-^^^' / dxu;i-x)u^{x) = 0. (10) 

J —00 

To study the consequences of ([T0|), we have to distinguish between three cases. 

If the two eigenvalues Ei and E2 are real, then E2 7^ Ei implies the vanishing of the 
integral on the left-hand side of (p!0[) , namely 

r+00 

dxul(—x)ui{x) = 0, (11) 

whereas for E2 = Ei = E, eq. ([ToD is automatically satisfied. In such a case, if E is 
nondegenerate, ^2(3;) is proportional to ui{x) = u{x) and following (H) we may propose 

/+00 
dxu*{—x)u{x) (12) 
-00 

as a counterpart of the standard normalization integral in PT-symmetric quantum mechan- 



ics. In the same way, the left-hand side of ([llD can be seen as a counterpart of the scalar 



product. It should be noted that in giving up the Hermiticity condition of the Hamiltonian 
and replacing it by PT symmetry, some properties of norms and scalar products have been 
lost. For instance, it is obvious that the integral in (|12]) is not positive-definite and may 
therefore be called pseudo-norm ||l7 . 



If instead only one of the eigenvalues (for instance £"1) is real and the other (£"2) is 
complex, eq. (|ll]) again follows from (plQl). 

Finally, if both eigenvalues Ei and E2 are complex, eq. (|T0|) implies condition (p!T| ) if 
E2 ^ E^. This is true in particular for E2 = Ei = E. In such a case, if E is nondegenerate, 
thence U2{x) is proportional to Ui{x) = u{x) and condition ([TT| ) becomes 

r+00 

dxu*{—x)u{x) = 0. (13) 



On the contrary, if E2 = = E* , then eq. (|T0|) is automatically satisfied and we may take 

/+00 
dx U2{—x)ui{x) (14) 

as a counterpart of the standard normalization integral. Note that in (0), ui{x) and U2{x) 
correspond to E and E*, respectively. It is remarkable that in going from real to complex 
eigenvalues, the roles of (|TTp and (|T^ are interchanged. Such results agree with the analysis 
carried out in [^. 

Let us consider in more detail the case of an eigenfunction u{x) corresponding to some 
real, nondegenerate eigenvalue E. It satisfies the time-independent Schrodinger equation 
— {h^/2m)(P'u{x)/(P'x + V{x)u{x) = Eu{x), and so is the case of u*{—x). Since E is nonde- 
generate, we must have u*{—x) = cu{x) for some complex constant c. From this relation, 
it follows that u*{x) = cu{—x) = \c\'^u*{x), hence |c| = 1. We conclude that 

u*{-x) = e''^u{x), < < 27r. (15) 

If = (resp. = vr) in (0), then u*{—x) = u{x) (resp. u*{—x) = —u{x)) or, in other 
words, u{x) is PT-symmetric (resp. PT-antisymmetric). If however 7^ 0, vr, it is always 
possible to convert u{x) into such a function by multiplying it by some appropriate phase 
factor. The functions 

t;<,(x) = eK<^-f +^^1 )m(x), a = ±l, (16) 

indeed satisfy the relations v*{—x) = av„{x), a = ±1 and are PT-symmetric (cr = +1) or 
PT-antisymmetric {a = —1). 

We propose to normalize u{x), or the corresponding Wo-(x), according to the prescription 

/+00 /■+00 
dxu*{—x)u{x) = a / dx[vi^{x)]'^ = a. (17) 
-00 J —00 

In the following, we illustrate the use of the above condition by considering some spe- 
cific cases of PT-symmetric model potentials, namely, PT-symmetric oscillator, generalized 
Poschl- Teller, and Scarf II potentials. In particular, we show that for this rule to be mean- 
ingful, we have to tell which PT-parity a, i.e., which function Vcr{x), we are going to associate 
to any given eigenfunction u{x) with real eigenvalue. 

5 



3 Modified Normalization for the PT- Symmetric Os- 
cillator 



The potential for the PT-symmetric oscillator is given by 

r>2 - i 

V^-\x) = {x-\cf + - r^, (18) 

[x — icy 

where a > 0. It can be obtained from the usual three-dimensional radial harmonic oscillator 
potential by effecting a complex shift of coordinate x — > x — ic, c > 0, and replacing / by 
a — |. Although it is beset with a centrifugal-like core, the shift of the singularity off the 
integration path makes the corresponding Hamiltonian exactly solvable on the entire real 
line for any a > 0, like the linear harmonic oscillator to which it reduces for a = \ and 
c = 0. 

If a is different from an integer (which we shall assume here), the PT-symmetric oscilla- 
tor Hamiltonian has a double series of energy eigenvalues ^ which may be distinguished 
by a quantum number q = ±1 For h = 2m = 1, they are given by 

= 4n + 2 - 2ga, n = 0,1,2, ... . (19) 

The corresponding eigenfunctions are expressible in terms of generalized Laguerre polyno- 
mials: 

u[l)(x) = ArJ^)e-5(-i^)'(x - ic)-''"+54-'?")[(a; _ ic)^]. (20) 

Here A^^^^) is some yet undetermined normalization coefficient. Although the Laguerre 
polynomials in (pOD are in principle defined only for —qa > —1, i.e., for a values in the 
interval < a < 1, their definition can be extended to any values of a. In the interval 
< a < 1 (and only in it), however, the energies ([l9|) are ordered according to increasing 
values of = 2n + |(1 — q) = 0, 1, 2, and alternatively correspond to g = +1 
and q = —1. In the special case of the linear harmonic oscillator (a = | and c = 0), q 
becomes the parity of the eigenfunctions. For this reason, q is called quasi-parity of the 
PT-symmetric oscillator. 
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It is straightforward to see that u^g^{—x) is related to as shown in (|T5D, where 

the phase is given by = n^—qa + |) — 2z/ (up to a muhiple of 2tt). Here u denotes the 
phase of the normahzation coefficient N'g"^ = \Afjjn^\e^'^. 

We now choose to identify the PT-parity a of the functions v„{x) in (0) with the 
quasi-parity q of the eigenfunctions (^). As a consequence, the functions 



gigf(-a+i) 



qn 



e-hi--^'^)' (x 



ic 



-qa- 



ic 



are PT-symmetric or PT- antisymmetric for g = +1 or g = — 1, respectively. 

Let us normalize the functions u'^'^{x) and according to the rule (|r 

is replaced by q. From (0), we obtain 



(21) 



where a 



dx 



,i7r(-i3Q!+^) 



•'^qn an i 



qn 



(22) 



where 



qn 



+ 00 



dxe 



-{x-icf 



IC 



-2qa+l 



{4-^")[(x-ic) 



(23) 



To calculate we consider the function of the complex variable z given by /(x 



qn 



. It is analytic everywhere except at z = 0, which is a branch 
point. Let us take the negative real axis as the branch cut. As the integral X^^-* can be 
seen as the limit for i? ^ oo of J^j^ dz f{z) for z = x — ic (i.e., the integral on segment AB 
in Fig. 1), we close a contour F = ABCDEF in the complex plane, avoiding the singularity 
by using a semi-circle Cp of radius p. Since f{z) is analytic within and on F, by Cauchy's 
theorem, we have Jj.dzflz) = 0. Evaluating this integral and taking the limits R —y oo 
and p — s> 0, we obtain that X^"^ converges if ga < 1, which is consistent with q = ±1 and 
< a < 1. 



As a result, eq. (|2^ ) becomes 

n+OO 



dx 



ATI") 



qn 



cos7r(— ga + ^ 



^ F(-ga + n + 1) 



(24) 



if < a < 1. Since the cosine on the right-hand side of (^4]) is positive or negative according 
to whether g = +lorg = — 1, the result in (|2^) is consistent with condition (|l^). The 
latter therefore leads to 



1/2 



As compared with the standard harmonic oscillator, the presence of a cosine on the right- 
hand side of (^) is a new feature. Such a cosine disappears for a = i, so eq. agrees 



with the known normalization coefficients of the linear oscillator (a = |, g = ±1) and of 
the three-dimensional radial oscillator (a = |, / = 0, g = —1). 

4 Modified Normalization for the PT- Symmetric 
Generalized Poschl- Teller Potential 



The PT-symmetric generalized Poschl- Teller potential is given by [|10 



V^(^-S) (a;) =[b^ + A{A + 1)] cosech^ t-B{2A + 1) cosech r coth r, r = x-ij, (26) 

where B > A + ^ > and — f < 7 < or < 7 < f . By using an sl(2,C) algebraic 
framework, it was recently shown [|1^] that whenever B — A — ^ is different from an integer, 
the corresponding Hamiltonian has two series of real energy eigenvalues Ej^^'^^ , which may 
be distinguished by a label q = ±1 (here called quasi-parity for convenience sake) and which 
for h = 2m = 1 are given by 



-S+n — ~ ~ 2 ~ ^) 1 ?1 — 0, 1, . . . , n+rnax, 

B-'l< n+^ax < 5 - I, (27) 

^(AB) ^ _ (^yi _ ^ n = 0, 1, . . . , n_jnax, 

A-l< n_^ax < A, (28) 
where B > ^ and A > 0. The accompanying eigenfunctions read 



= Arj^'^)(y-l)^(^+i)(y+l)i(^+|)pf'/^)(y) 



7/ 

(jn V"^y " " qn 



= A/JAB)2^(A+A«+i) (sinhi) (cosh l)''^ P(^'^) (cosh r), (29) 

where A = q{A — S + |), /i = —A — S — |, ?/ = coshr, P^^'^\y) is a Jacobi polynomial, 
and Mjjn'^^ some normalization coefficient. 

In the limit 7^0, the real generalized Poschl- Teller potential is recovered and only 
the eigenvalues E^^^ and corresponding eigenfunctions m^'^^(x) survive Since this 



potential is singular, it has to be restricted to the half-line (0, +cxd). The complexified 
potential, as given above, gets regularized on performing the shift x — x — 17 and so may 
be considered on the entire real line, as it was the case for the PT-symmetric harmonic 
oscillator in the previous section. 

The analogy between these two PT-symmetric potentials also extends to the normaliza- 
tion problem. The phase in (|T5|) is now given (up to a multiple of 2tt) by = 7r(A + 1) — 2z/, 
where u denotes the phase of the normalization coefficient. This shows that the role of a for 
the PT-symmetric oscillator is now played by the constant —qX = B — A — ^. As in sect. 3, 
the PT-parity a may be identified with q, leading to PT-symmetric and PT-antisymmetric 



eigenf unctions v^^'^\x) (g = ±1), given by 



j^(A,B)^ i^ginh f )^^^ (cosh lY^^' ''^^(coshr). (30) 



Applying again prescription (0) to normalize the eigenf unctions, we are led to the integral 

P^^''^)(cosh z) on the contour P of Fig. 1, 



z\2\+l 

2' 



(coshD^A'+i 



of the function f{z) = (sinh 
where c is replaced by 7 (and the latter is assumed positive for simplicity's sake). In this 
respect, it is worth noting that among the branch points of f{z) at 2; = 0, iivr, ±2i7r, 
±3i7r, . . . , we only have to take care of the first one because of the range of 7, as stated 
below eq. (|26|). 

The final result for the normalization integral reads 



dx 



u 



(A,B) 
qn 



i-x) 



qn 



COS 



vr(A + 1)1 /(^'^^ 



(31) 



provided A is restricted by the condition A > —1, which is consistent with q = ±1 and 
A + I < B < A+ I. In (|31]), 4^''') is the real integral 



/{a,m) 



dt {t-i)\t + iY \Pi^'''\t) 



(32) 



which, for g = — 1, appears in the calculation of the normalization coefficient of the real 
potential eigenfunctions. It is convergent and positive for the above-mentioned restricted 
range of parameters and both values of q. As far as we know, no closed analytical formula 
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is known for it, except when n = 0. In such a case, the change of variable s = 1/t leads 



to 23 



j(A,/.) ^ oA+M+i r(-A-/i-i)r(A + i) 
° r(-/i) 

Furthermore, the cosine on the right-hand side of ( PT] ) is positive or negative according to 
whether g = +lorg = — 1. Hence eq. (^) is consistent with from which it follows 
that 

Ar(„^-^)| = {2cos[7r(A - 5 + 1)]/^'^)}"'^' . (34) 

The cosine disappears from ii B = A + 1. For this choice of parameters, potential (P^) 
reduces to V'''^'^^^^ = —^{A + l)(2y4 + 1) sech^ |, which remains nonsingular in the limit 
7^0. 

As a final point, it is worth mentioning that the generalized Poschl- Teller potential (^) 
is related to the PT-symmetric Poschl- Teller II potential considered in |Q through a com- 
plex point canonical coordinate transformation [^, which also connects the corresponding 
eigenf unctions. The results of the present section can therefore easily be applied to the 
latter potential. 

5 Modified Normalization for the PT- Symmetric 
Scarf II Potential 

The PT-symmetric Scarf II potential is defined by [p!0[1 

V^^'^\x) = - [b^ + A{A + 1)] sech2x + i5(2A + l)sechxtanha;, (35) 

where A > B — ^ > 0. It was shown using sl(2,C) as a tool [|lOl that for values of 
A — B + ^ different from an integer, the corresponding Hamiltonian has a double series of 
energy eigenvalues Ejj^'^\ distinguished by a label g = ±1 (again called quasi-parity for 
convenience sake), and given by 

^11'^^ = -{A-nf , n = 0, 1, . . . , n+^ax, 

A-1 < n+rnax < A, (36) 
-Sil''^^ = -[B -\-ri^ , n = 0, 1, . . . , ri-rnaxj 

B-'l< n_^ax <B-\, (37) 
10 



a h = 2m = 1, A > 0, and B > \. The accompanying eigenfunctions read 

= A/;(^'«)(sechx)-5(^+^+i)e"^(^-^)"'^^*'^"(""^")pf '^^ism^ (38) 

in terms of Jacobi polynomials. Here A = q[—A + i? — |), /i = —A — i? — |, and M^^'^^ is 
some normalization coefficient. 

Contrary to the potentials considered in the two previous sections, the PT-symmetric 
Scarf 11 potential is not obtained from the corresponding real potential |^ by a com- 
plex shift of coordinate, but instead by complexifying one parameter [B iB). Such a 
complexification is responsible for the appearance of the additional series (|37D of energy 
eigenvalues. 

Comparing ( |38D with (|15|), we find (up to a multiple of 27r) = — 2z/, where v is the 
phase of the normalization coefficient. If we identify PT-parity with quasi-parity again, we 
obtain 



M^^'^^ (Sech X)-^ (^+^+^) e" ^ (^-'^) arctan(sinh X) p^A,M) ^ j gj^j^ ^ ^gg^ 



where exp[i(g — 1)^] = 1 for PT-symmetric functions and — i for PT-antisymmetric ones. 
The normalization integral is now 



dx 



qn 



dx (sech X 



|— A— /i— Ig— i(A— ^) arctan(sinhx) 



P^^'^)(isinhx) 



(40) 



Let us consider, for instance, the n = case corresponding to the lowest states of both 
series of energy levels. The integral on the right-hand side of (^) is then easily calculated 
by performing the change of variable tany = sinhx. The results read p3[ 



+ 00 

-oo 

+ 00 

-oo 



dx 
dx 



U 



(A,B), 
+0 I 



-X 



U 



+0 



X] 



Mir' 



TxT{2A) 



2^A-^V{A-B+\)T{A + B + \ 
t:V{2B - 1) 



TY' (41) 



-V (42) 



2^B-2Y(^B-A-\)T{B + A+^, 

The right-hand side of (PD is positive for any allowed values of A and B because 
the arguments of all three gamma functions are positive. This is compatible with the 
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normalization condition (|T^ with a = q = +1. On the right-hand side of (^), on the 
contrary, the argument of T{B — A — ^) is negative while those of the two remaining 
gamma functions are positive, which means that the sign of the normalization integral for 
(7 = q = —1 depends on the relative magnitudes of A and B — ^. A negative sign compatible 

with (|1^ is obtained ioi B - ^ + 2k < A < B + ^ + 2k, k = 0, 1, 2, In particular, for 

= 0, we get a condition B — ^<A<B + ^ rather similar to the conditions previously 
encountered for the PT-symmetric oscillator and generalized Poschl- Teller potentials. 

6 Conclusion 

To conclude, we have derived, in the set up of a very general framework, the conserva- 
tion law pertaining to PT-symmetric quantum mechanical systems. We have found, as a 
consequence, that the normalization integral ceases to be always positive definite. The pos- 
sibilty of encountering non-positive definite norms is illustrated by means of some specific 
examples of PT-symmetric potentials. 
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Figure caption 

Fig. 1. Contour F in the complex plane. 
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